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Abstract 

We present a full introduction to the recent devised perturbation 
theory for strong coupling in quantum mechanics. In order to put 
the theory in a proper historical perspective, the approach devised in 
quantum field theory is rapidly presented, showing how it implies a 
kind of duality in perturbation theory, from the start. The approach 
of renormalization group in perturbation theory is then presented. 
This method permits to resum secularities in perturbation theory and 
makes fully algorithmical the resummation, transforming the pertur- 
bation calculations in a step by step computational procedure. The 
general theorem on which is founded a proper application of the strong 
coupling expansion, based on a result in the quantum adiabatic the- 
ory, is then exposed. This theorem gives the leading order of a strong 
coupling expansion. Then, after the introduction of the principle of 
duality in perturbation theory that puts in a proper context the quan- 
tum field theory method, the resulting theory of the strong coupling 
expansion and the free picture are presented. An algorithm for the 
computation of the perturbation series is finally given. This approach 
has a lot of applications in fields as quantum optics, condensed matter 
and so on, extending the original expectations of the quantum field 
theory method. So, we give some examples of application for a class 
of two-level systems that, in recent years, proved to be extremely im- 
portant. One of the most interesting concepts that can be obtained in 
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this way is that of a Quantum Amplifier (QAMP) that permits to ob- 
tain an amplification to the classical level of the quantum fluctuations 
of the ground state of a single radiation mode. 
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1 Introduction 



After the discovery of the quantum chromodynamics (QCD) and asymp- 
totic freedom 1 , it become increasingly important to treat in some 
way a theory having a non-perturbative behavior. Since then, different 
approaches have been devised to recover the spectrum of QCD at low 
energies but here we focus on a perturbative method that, although 
did not prove to be useful to treat QCD problems, paved the way 
toward a strong coupling expansion with a possible wider scope. This 
approach, the strong coupling expansion, reached its best formulation 
in a paper by Bender and coworkers El El where, being applied to 
a X(ft 4 quantum field theory, it was proved that a lattice formulation 
could give a manageable formulation but that the problem is moved on 
taking the limit of zero lattice spacing on some very singular series for 
the relevant quantities of the theory. Some resummation techniques 
were devised without much success. 

In quantum field theory it is customary to start with a path integral 
in the Euclidean space like for a X(ft 4 theory (here and in the following 
h = c = 1)0 
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and to obtain a weak coupling expansion in A as 
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where we have extract formally the quartic term from the path integral 
and put a functional derivative for each power of the field. Then, we 
are left with a gaussian integral that can be easily computed and we 
have the sought expansion that can be cast in the form 



oo 
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At this stage one may ask what happens if we do the opposite oper- 
ation, that is, if we consider the kinetic term ^(d(ft) 2 + ^m 2 (ft 2 as e 
perturbation. This can be realized rewriting the path integral as 
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being G~ 1 (x—y) = (— d 2 +m?)5{x— y) the inverse of the free Euclidean 
Green function. It can be proven that, on a lattice being the theory 
highly singular, the expansion takes the form 2 



where one can see a dependence on the inverse of the coupling con- 
stant A. We see that we have arbitrarily chosen different part of the 
Euclidean action as a perturbation, a possibility offered by the free- 
dom proper to this choice, and, by doing that, we get two perturba- 
tion series having as development parameter one the inverse of the 
other. This kind of "duality" can be found e.g. in fluid mechanics 
with the Navier-Stokes equation where we can have different pertur- 
bative regimes by taking as unperturbed part the Eulerian term or 
the Navier-Stokes term 5 . This regimes are characterized by large 
Reynolds number and small Reynolds number respectively. This kind 
of duality in perturbation theory can be indeed seen rather widely 
and appear an ubiquitous property due to the freedom in the choice 
of what a perturbation is. The first formulation in this sense in quan- 
tum mechanics appeared in Ref.jO] following a series of works where a 
strong perturbation theory in quantum mechanics has been formulated 
QISI03 E3 • In this work it was shown that the adiabatic approxima- 
tion, formulated as in [111 I12j . is the leading order approximation of 
what can be called a dual Dyson series, as the weak coupling expansion 
is the well-known Dyson series. 

Being a dual series to the standard Dyson expansion it shares the 
same problems. Particularly, one of the most relevant questions one 
has to face is that of secularities ^JJEJEIl- ^ secularity is a polyno- 
mial contributions to the series that increases without bound making 
the series itself useless. The name "secularity" is taken from celestial 
mechanics where firstly these terms appeared in perturbation series 
with a secular timescale. Then, in order to have an useful tool to 
make computations one has to devise a way to remove such singular 
terms. In the course of time several approaches have been proposed to 
this aim |13[ I14 | IT5]. The difficulties with these methods are essentially 
linked to the impossibility to make them algorithmic in some way. But, 
recently a new approach has been proposed ^f3ElE2E3 that solved 
this problem making perturbation computations straightforward to re- 
alize. This method relies on the renormalization group techniques that 
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aim to find the envelope of the computed perturbation series to the de- 
sired order. This approach has been successfully applied in quantum 
mechanics and to the dual Dyson series [!M ETT1 1221 [21 [25] . 

The existence of a strong coupling expansion can prove to be very 
important as it gives the opportunity to study the solution of a differ- 
ential equation in different regimes in the parameter space. In turn, 
this means that new physics can be uncovered. Although the strong 
coupling expansion has found several applications in different fields as 
strong atom-laser interaction [^Hl H7j , quantum chaos |2H1 HH1 EH] and 
quantum Zeno effect [HD 02] , the workhorse for a lot of studies in fields 
as quantum optics, quantum computation and condensed matter is the 
two-level system j^HUM]. Then, to present a lot of examples of appli- 
cation of the strong coupling expansion, we analyze some examples of 
two-level Hamiltonians showing how relevant physics can be exploited 
in this regime giving finally rise to the concept of a quantum ampli- 
fier (QAMP), able to amplify the quantum fluctuations of a radiation 
mode in the ground state to the classical level jHSl HUH E3 ESI • 

This review is structured as follow. In Sec|21 we introduce the 
duality principle in perturbation theory for the Schrodinger equation. 
In SecEJthe fundamental theorem of the strong coupling expansion to 
the leading order is proved. In Sec|I]the problem of secularities is put 
forward as a general problem in perturbation theory that is inherited 
by the strong coupling expansion. In SecEJthe renormalization group 
method to solve the secularity problem is presented. In SecEl the 
approaches discussed above are merged to formulate an algorithm for 
doing perturbation theory. In Sec|7| we apply the strong coupling 
expansion to a two-level model well-known in quantum optics and 
finally, in sec|Sl we show how the two-level model discussed in Seed 
can amplify the quantum fluctuation to the classical level producing 
a classical field (QAMP). In SecElthe conclusions are given. 

2 Duality in Perturbation Theory 

The starting point of our analysis is given by a quantum system with 
a Hamiltonian 

H = H + V. (6) 

It is usually assumed that the dynamics for th Hq part is known. 
Adding the V part can make the problem unmanageable unless one 
uses perturbation theory and the V part is smaller in some sense 
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with respect to Hq. To account for this, we introduce an arbitrary 
parameter A that we now consider small but we take to be unity at 
the end of the computation. So, we write H = Hq + XV and aim to 
solve the Schrodinger equation 



(H + \V)\iP) =i 



dt 



(7) 



When A is a small one does the transformation (interaction picture) 



obtaining the equation to solve 



-M t 



dm 

dt ' 



(8) 
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The solution of this equation, generally know as Dyson series, can be 
written as 



\ipi) = Texp 



-i [ dt't 
Jo 



Mot' 



We 



-iH t' 
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being T the so called time ordering operator and \ip(0)) the initial wave 
function. This is just a formal writing for the small perturbation series 



M)=\I- 
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dt'e 
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Jo 



i I dt'e 1 
'o Jo 

(11) 

that we recognize as a power series, \ip[) = J2^=o ^ n |^n(i))) that can 
have a meaning only for small values of A. 

We have done some assumptions to derive the Dyson series. We 
have assumed that Hq is the unperturbed part of the Hamiltonian 
whose dynamics is known and that A was a small parameter. Then, 
we have obtained a result that could be meaningful, at least asymp- 
totically. But we can relax both these assumptions in view of the fact 
that the choice of an unperturbed part and a perturbation is totally 
arbitrary and one may ask what would be the face of a series where 
the role of Hq and V is interchanged. We can easily work out this 
exchange into the Dyson series obtaining 



We 



-Mot" 



+ 



\4) F ) = Texp 



: I dt' 

Jo 



JXVt' tt -i\Vf 
e tiQe 



W0)>. 



(12) 
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To understand what we have done we also transform the integration 
variable as r = Xt' and we have 



that means the series 



1 
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X 



Xt 



dT'e iVr 'H e~ iVT ' + • 



(14) 

A meaning can be attached to this series only if we take the limit 
A — ► oo and we have obtained a dual series with respect to the Dyson 
series having a development parameter inverse to the latter. So, using 
a symmetry in the choice of what a perturbation is we were able 
to uncover a new perturbation series, dual to the small perturbation 
series, useful in problems with a strong perturbation. Finally, we can 
set A = 1 and we will have the Dyson series as 



\ipi) = Texp 



: f dt't 
Jo 



,iH t' 



Ve 



iH t 



(15) 



after passing to the interaction picture with the unitary transforma- 
tion ex.p(—iHot) and a dual Dyson series 



■ fdt'j vt ' 
Jo 



H Q e 



-iVt' 



Mo)). 



(16) 



passing to the free picture with the unitary transformation exp(— iVt). 
It is easily realized that the two series will share the same kind of 
problems as secularities, asymptotic convergence or divergence and so 
on. But, we are aware that this does not diminish the usefulness of 
the perturbative approach. 



3 The Leading Order of the Strong 
Coupling Perturbation Theory 

The leading order of the dual Dyson series is quite straightforward 
to define if the perturbation V does not depends on time. Things 
are more involved otherwise. Indeed, the dual Dyson series must be 
redefined if V depends on time and we will prove that a series, derived 
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formally from the adiabatic approximation, due to Mostafazadeh 
is recovered in this case. Again we consider the problem 

(H + XV(t))\m)=i^^- (17) 

with A — > oo. To recover the dual Dyson series, we have to determine 
the unitary operatore Up(t) such that 

\V(t)U F (t) = i^^ (18) 

with the same limit for A. It is very easy to recognize here the starting 
point of the proof of the adiabatic theorem in quantum mechanics 38 , 
12 and we are recovering a series, in the same framework, due to 
Mostafazadeh |11| I12|. This is due to the slowing down implied by 
the parameter A going to infinity. Then, without recurring to any 
adiabatic hypothesis, we can write [""] 

U F (t) = Y J ^ n(t) e- iI o Vn(t ' )dt '\n;t)(n;0\ (19) 

n 

being 7 n (t) the geometric part of the phase and v n (t) the dynamical 
part such that V(t)\n;t) = v n (t)\n;t). We will show that this gives 
the proper dual Dyson series. Higher order corrections can be written 
down and are given in Refs. [3*51 fTTl "T2*] . It is interesting to note that 
this theorem has found an application in the studies of Zeno effect in 
quantum systems [31] . 

Once the unitary operator Up(t) is known, we are able to pursue 
the computation to the end. Indeed, we have to solve the Schrodinger 
equation 

H F \Mt))=i 9 -^ (20) 

being 

Hp = ^^e i [ 7 "W- 7m Wl e - i J'o[ 1 '"W- 1 ' m W](m,t|^o|n,t)|m,0)(n,0| 

m n 

(21) 

the transformed Hamiltonian. We easily realize that the Hamiltonian 
Hp can be split in two parts as 

Hp = ^2(n,t\H \n,t)\n,0)(n,0\ (22) 

n 

+ ^ e i ^W^ m ^e^^^^-'^(*')] <e/ (m,*|fl'o|n,t)|m J 0)(n,0| 
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and we are able to obtain the analogous equations of the interaction 
picture for the probability amplitudes that apply in the case of the 
dual Dyson series. In order to obtain this result we write the solution 
of ea.pU]) as 

\M*)) =J2 C n( t >~ if ° h ° n{t ' )dt '\ n ^ ( 23 ) 
n 

being 

h 0n {t) = (n,t\H \n,t). (24) 
This gives the equations for the amplitudes 

ic n (t) = ^e < [' 5 *W-'> fe W]e- < -'"o[ e k(<')- e »«')]*'(n 1 t|flo|fc,t>Cfc(t) (25) 
k 

with e n (t) = ho n (t) + v n (t) . These equations are similar to the ampli- 
tude equations in interaction picture, normally found in textbooks. 

Some considerations are in order at this point. By the adiabatic 
theorem we can only apply the above approach for a discrete spectrum 
on the perturbation. This in turn implies that, if the spectrum of the 
perturbation is continuous as in the coordinate space, then a lattice 
regularization is needed, This takes us back to Bender et al. approach 
PI El E] with all the difficulties this means by taking the limit of the 
lattice spacing going to zero. 

4 The Secularity Problem in Pertur- 
bation Theory 

A perturbation series is plagued by secularities when polynomial terms 
in time appears in it. These terms have the property of being not 
bounded for large time making the series generally useless. 

In order to have an idea of what really happens in these situations, 
let us consider the well known quantum mechanical problem 

H = -^-03 + 9&i cos(cjt) (26) 

representing a two level atom driven by an oscillating field with fre- 
quency uj. 01,03 are Pauli matrices, g is the coupling constant and 
A is the separation between the atom levels. This problem has a 
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large body of literature due to its vast field of applications (see e.g. 
Ref.JUllISl)- 

We can apply to this problem all the machinery devised in the 
preceding section for the strong coupling expansion. So, one has 

U F (t) = e -^ s »M)| + )(+| +e i£sin(«,t)|_^_| (27) 

being ci|±) = ±|±). By using the explicit expression for the states 
|±) one can prove that U F (t) = e~ ieT1 i sin(arf) as it should be for the 
free picture. It is very easy to obtain 



H F= * 
2 



,-<^sin(wt)|_^ + | +e *^sin(wt)| + ^_|j _ ^8) 



Looking for a solution in the form \ipp(t)} = c+(t)\+) + c_(t)|— ) we 
obtain the equations for the amplitudes 

ic+(t) = ^e^^c-^) (29) 

id_(t) = ^ e -^ sin M) c+ (t) 

that give rise to the perturbation series till first order 

c+(t) = c + (0)-i-J o (z)c-(p)t--Y,M*) c_(0) + (30) 

2 2 *—f nu 

A A * v e —inu>t _ j 

C_(t) = C_(0)-Z-J O (*)c + (0)t + - V Jn(z) C+(0) + --- 

2 2 *—f nw 

where use has been made of the relation e tzsm ( ujt ) = J n (z)e mult 
with J n (z) the n-th Bessel function of integer order and z = ^ in our 
case. We see immediately that the perturbation series is plagued with 
secularities and so is useless at this stage. We have to understand 
what is going on here by properly resum such terms. In this way, 
we will discover here a physical effect, i.e. Rabi oscillations between 
the states |±). The resummation technique to accomplish our task is 
described in the next section where we will complete our computation. 



5 Renormalization Group Method for 
the Resummation of Secular Terms 

The method of renormalization group to resum secularities in a per- 
turbation series was firstly proposed in |16U17| . Here we present an el- 
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egant reformulation obtained by the mathematical theory of envelopes 
by Kunihiro [T31IT3]. 

Kunihiro approach can be described as follows. Let us consider 
the following equation 

x(t) = f(x(t),t) (31) 

being x(t) a vector in M n . The initial condition is given by x(to) = 
X(to). At this stage we assume X(to) not yet specified. We write the 
solution of this equation as x(t; to, X(to)) that is exact. If we change 
to to t we are able to determine X(t ) by assuming that the solution 
should not change 

x(t;t ,X(t ))=x(t;t' ,X(t' )) (32) 
that in the limit to — * t'o becomes 

dx _ dx dx dX 

dt dt dx dt { ' 

giving the evolution equation or flow equation of the initial value 
X(to). We recognize here a renormalization group equation and this 
gives the name to the method. 

Till now, all our equations are exact and no perturbation theory 
entered in any part of our argument. But, except for a few cases, 
the solution x(t;to,X(to)) is only known perturbatively and such a 
solution are generally valid only locally, i.e. for t ~ to and t ~ t' and 
a more restrictive request should be demanded to our renormalization 
group equation 



dx 
dt 



dx 

t =t dt o 



dx dX 

t0=t + dxW 



= 0. (34) 

t =t 



But this equation can be interpreted by the mathematical theory of 
envelopes |18| . Indeed, varying to we have that x(t; to, X(to)) is a fam- 
ily of curves with to being a characterizing parameter. Then, eq.(J3U) 
becomes an equation to compute the envelope of such a family of 
curves. Such an envelope is given by x(t;to = t) = X(t), the initial 
condition. It can be proven that X(t) satisfies the equation (|31[) in a 
global domain up to the order with which x(t; to) satisfies it locally for 
t ~ to- This gives the condition for the computation of the envelope 



,lr = (35) 
t =t 



dt 
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The Kunihiro method is very effective to build resummed pertur- 
bation series, eliminating the seculairities that appear to plague them. 
Besides, it permits to transform a perturbation computation in an al- 
gorithm straightforward to apply as we are going to see in the next 
section. 



6 An Algorithm for Doing Perturba- 
tion Theory in Quantum Mechanics 

In order to exploit what we mean by an algorithmic computation of a 
perturbation series, we come back to the example given in sec|IJ All 
we have to do, as our first step, is to recompute the perturbation series 
at a generic initial time to and to assume generic initial conditions. 
This yields 

c+(t) = c + (t )-i-J (z)c-(t )(t-t )--J2Jn{z) c_(t ) + --{36) 

2 2 nuj 

c_(i) = c^(t )-i-Jo(z)c + (t ){t-t ) + -Y^J n {z) £+.(to) H ■ 

2 2 z — ' nuj 

We realize easily that the envelope could not be computed with this 
series as is. What we need here is to dress all the phases in the 
exponentials. This gives 

^ ^ e inuit _ e -inu)4>(to) 

c+(t) = c+{t )-i-J (z)c-(t )(t-t )--y2Jn(z) c_(t )+(37) 

2 2 ^-^ nuj 

^ A e ~inu)t _ e inuj4>(to) 

c_(t) = c-(t )-i—Jo(z)c+(to)(t-to) + — VJ„(z) c+(t ) + ---. 

2 2 ntu 

where we have introduced a renormalizable phase </>(to) = —to- This is 
strictly linked to the property of quantum systems to have a freedom in 
the choice of the initial phase. At this point we use the renormalization 
group or envelope equation (|33)) giving 

^± + l ^J ( z )c_(t) + ... = (38) 

dc-(t) A T . . . , 

' +i-Jo{z)c+{t) + ... = 



dt 2 

mi 

dt 
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+ ... = 



and the perturbation solution is then given by 



A e inujt _ 2 

c+(t) = £ f ( i )-_£j n (s)___eL( f ) + ... (39) 

A , e -inwt _ j 
C_(i) = c_(t) + -2j JnW C+(t) + ---. 

completely solving our problem, till first order, and having all secular 
terms properly removed in an algorithmic and simple way. 

We can finally exploit fully our algorithm for doing perturbation 
theory in presence of secular terms in the strong coupling regime for 
quantum mechanics. The rules for computing the unitary evolution 
operator, at any desired order, are the following |23| : 

1. Consider the following unitary transformation on the equation 
U F (t) = ^e i ^ {t) e- i ti Vn{t ' )dt '\n-t)(n-U\ (40) 



with the eigenstates of the perturbation \n;t). This gives the 
transformed Hamiltonian 

H F (t) = U F (t)H U F (t). (41) 

The dual Dyson series is computed by [3 |H1 El E3 



5d(Mo) = ^exp 



-te 



H F (t')dt' 



t 



(42) 



being as usual T the time ordering operator and an ordering 
parameter e has been introduced that will be taken unity at the 
end of computation. It is fundamental for our argument that 
the computation of this series is performed at a different starting 
point to- 

2. Assume, at the start, that the time evolution operator has the 
form 

U(t, t ) = U F (t)S D (t, t )U R (t ) (43) 
where Un(to) is a "renormalizable" part of the unitary evolution. 
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3. At the given order one gets Spit, to) as 

S D {t, t ) = I - iefi(t, t ) - e 2 f 2 (t, t ) + ... (44) 

and, at this stage, if some oscillating functions in to appear like 
e~' lujt ° then introduce the phase <j)(to) = —to as a "renormaliz- 
able" parameter rewriting it as e^i 1 ®) . The secularities must be 
left untouched. 

4. Eliminate the dependence on to by reauiring jlHl I19j 

dU(t,t ) 



= (45) 

t =t 



dt 

and one obtains the renormalization group equations 

= e 9l U R (t)+e 2 g 2 U R (t) + 0(e 3 ) (46) 



dt 

e(/>i^(t) + e^ 2 0(t)+O(e J 



d ^ -± ±r^ , 2 , , U \ , rM 3^ 



dt 

where, at some stage, to obtain such equations at the second 
order, we have to use their expressions at the first order as, to 
compute their form at order n-th one have to use these equations 
at the order (n — l)-th, into the condition (|45|). This is a step 
toward the computation of the envelope of the perturbation series 
as said in secEI 

5. Finally, the renormalization equations should be solved and sub- 
stituted into the equation 

U(t,t )\ t0=t (47) 

giving the solution, i.e. the envelope, we were looking for without 
secularities at the order we made the computation. 

Once the unitary evolution is known, we can easily compute the 
wave function, given the initial condition, recovering the case we have 
shown of the driven two-level system. 

We are going to see this approach at work in the next sections. 

7 Two-Level Systems and the Strong 
Coupling Perturbation Theory 

In quantum optics the interaction between a single radiation mode 
and a two-level atom proves to be a paradigm for most of applications 
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[33 HaE]. The Hamiltonian is 

A + + 

H= —a z +uja) a + ga x {a^ + a) (48) 

that differs from Hamiltonian (|26|) by having a fully quantized radia- 
tion field of frequency u>, rather than a classical field, whose creation 
and annihilation operator are a* and a. 

The standard approach ""[Bj to this problem is given by doing the 
unitary transformation 

and then the rotating wave approximation is applied keeping only the 
near resonant terms (A ~ to). This reduces our model to the well- 
known Jaynes-Cummings Hamiltonian |39l I4U| 

Hjc = -j°z + g(acr + + a) a J) (50) 

being now 5 the detuning between the frequency of the radiation field 
<jj and the separation between the levels of the atom A. This aspect 
is well known having a large body of literature since its inception and 
being a foundational matter for quantum optics. 

Instead, here our aim is to realize a complete study of the Hamil- 
tonian Q48|) from the point of view of our strong coupling approach 
[121 14*3*1 EU EH] • F° r our aims we have to compute U p (t) and this is 
done by solving the eigenvalues problem 

[loo) a + ga x (a) + a)]\[n; A]) = E n>x \[n; A]) (51) 

whose solution is given by 

|[n;A])=e^- at )|n)|A) (52) 

with a x \X) = A|A), A = ±1 and the eigenvalues, independent on A 

2 

being E n = nu — Then, it is straightforward to write the unitary 
evolution operator as 

U F0 (t) = ^2e- iE - t \[n;X])([n;X]\\X)(X\ (53) 

that gives rise to the Hamiltonian 

Hf = U F0 (t)-a z U F0 (t) (54) 
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It is easily realized that it can be rewritten in the form |42l 1431 IHl I45j 

Hp = Hq + H x (55) 

being 

K = jT, e ~^ L - (S) [IM]><[™;-i]l|i)<-i| + l[n;-i]><[n; i]||-i>< 

(56) 

being L n the n-th Laguerre polynomial [IS] an d 

(n\e-^ a -^\m)\[n-,l])([rn-,-l]\\l)(-Xm 



H 



-i(n—m)ujt 



(ri|e^^- at )|m}|[n; -l]){[m; 1]|| - 1)(1| 



At this point we can iterate the procedure by diagonalizing the Hamil- 
tonian Q56|) . The eigenstates are 



|^;(7> = -^H[n;l]>|l> + |[n;-l]>|-l>] 



and the eigenvalues are 



~ A -2sl T (kg 2 



being a = ±1. So, we can write the unitary transformation 
U F1 (t) = Y,e- l ^ t \^n;a}(iP n ;a\ 



11.(7 



and get the transformed Hamiltonian 

H{ = ul^HxUpxit) 

that is 



(58) 



(59) 



(60) 



(61) 



H 'l = T E E ^n,, 1(r2 e- i[(,l - m) ^ (l "^ 1 -^ V2)1 Vn^l>(^m;<7 3 | 



m,n,m^n o"i,<T2 



being 



l r 



m)£j 2 



(62) 
(63) 
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So, we have accomplished the unitary transformation 

U F (t) = U F0 {t)U Fl {t) (64) 
and we are left with the Schrodinger equation for our aims 

HiS D (tM=i^^ (65) 

that we solve by the perturbation theory obtaining the strong coupling 
expansion for this problem. As said in the formulation of the algorithm 
in Sec® we assume a solution in the form 

U(t, t ) = U F (t)S D (t, to)U R (to) (66) 

being Un(to) a renormalizable part of the unitary evolution. 

In order to obtain the sought series we need to understand where 
resonances occur, that is where the condition (n — m)uj — {E n ^ ai — 
E m .a 2 ) = is met. This happens for n ^ m and for the two other 
conditions o\ = 02 (intraband resonance) or o\ ^ 02 (interband res- 
onance). Then, as required by our algorithm, we compute the dual 
Dyson series at an initial time to obtaining 



SD(t,t ) = I ~ i ~2 



\ipn]cr 1 )(Tp m ;cr2\(t-to) (67) 



intraband 



+ fcmn,aia 2 \' l Pn,(Tl}{' l Pm)<T2\(t-to) 



interband 



m,7l,CTj ,(T2 

out of resonance 

rn)u>-(E„, CT1 — £ m , CT2 )]i _ -i[( y n-m)u]-(E n ^ 1 -E m ^ 2 )}t 



-i[{n - m)uj - (E n>(T1 - E m ^)\ 

The next step is to introduce the phase 4>{to) = —to into the exponen- 
tials changing the series into 



-|V>n;<7i)(V> m ;o"2| + 



\tpn;^i){tp m ;cT2\(t - 1 ) 

Proband 

+ ^ ^■mn,a 1 cr 2 IV'n;0-l>(V'm;0-2|(i - ^o) 



interband 



m,n,<T j ,o"2 

out of resonance 
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e —i[(n—m)ui—(E n , C F 1 —E mtC , 2 )]t _ e i[(n-m)u)-{E n ^ 1 -E m ,a 2 )\4>{ t o) 

-i[(n - m)ijj - (E n ^ - E m ^ 2 )] 

Finally, we can compute the envelope of U(t,to) obtaining the renor- 
malization group equations 



mi 1 



+ 



dU R {t) 
dt 



d<j){t) 
dt 



-i — 
2 



mi 0"2 



.intraband 



interband 



mi f 2 



+ ... = 0. 



Then, computing U(t,to)\t =t, we get the series 



U(t) 



U F (t) 



1 + 



A 



out o/ resonance 



e -i[{n-m)u)-{E n ^ 1 -E m ^ 2 j\t _ ^ 

(n - m)w - (-En.CTi - E m . a2 ) 
U R (t). 



(69) 



t/ fl (t) + . . . 



(70) 



The relevant result is that, by the renormalization group method, 
we have resummed the perturbation series obtaining the unitary evo- 
lution, Un{t), proper to Rabi oscillations as it should be 
plus a first order correction. Rabi oscillations in the strong coupling 
regime, as described here, have been recently observed in Josephson 
junctions .47] • 



8 An Application: The Quantum Am- 
plifier (QAMP) 

The next step is to generalize the model of sec|7]to N two-level atoms. 
We will find a new physical effect that can be seen as a quantum 
amplification of the vacuum fluctuations, that is, we realize a quantum 
amplifier (QAMP). For our aims, it is very easy to generalize the 
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Hamiltonian ()48j) as 

N N 

2 



Hn = —^2cr Z i + u}a i a + g^2a x i(a ! + a). (71) 

i=l i=l 



Now, we introduce the analogous of angular momentum operators as 

1 * 

i=l 

1 * 
= 2 E » 

8=1 

1 * 

^ = TjE ^ 

i=i 

q2 q2 i q2 I q2 

D — D x -\- Oy -\- D z 

with the well-known commutation relations [S, Si] = and [Si,SA = 
ieijkSk, with the index k that can take the values x, y, z. Now, we 
have that, depending on N being even we can have a zero momentum 
state, otherwise [IB] 

\S X \ <= S <= y (73) 
N N 

-y <= s ' <= t 

giving the Dicke states \S, S x ). 

At this stage, we can iterate the procedure in sec|7]by diagonalizing 
the Hamiltonian 

H F = uja)a + 2gS x {a) + a) (74) 

with the eigenstates 

\[n;S,S x }}=e^ a - a ^\n}\S,S x } (75) 

and eigenvalues 



n - - 9 X 



uj (76) 



and this time we have no degeneracy with respect to the Dicke states 
as happened for the single two-level atom. Then, it is straightforward 
to write down the unitary transformation as 

U Ht) = E E e- iEn ' sJ \[n; S, S x ])([n; S, S x ]\. (77) 

n S,S X 
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Already at this stage we can have quantum amplification. Indeed, 
let us take as initial state \ip(0)) = ^hat * s > we nave ^he 

radiation field in the ground state and the maximal Dicke state. In 
the ground state, the radiation field has vacuum fluctuations as it is 
well-known. This gives 

n 2 2 Vn! 

that is nothing else that the solution of Ref. [351 071 13*3*] . that is, a 
coherent state with a parameter increasing for large N. At this stage, 
we can take two different thermodynamic limits. The first one is given 
by statistical mechanics, i.e. ./V — > oo, V — > oo and y = const, being 
V the volume that contains the radiation mode (e.g. a cavity). The 
second one is given just by the limit N —* oo keeping fixed the volume. 
In the former case, we observe that g cx -if- and so, the thermodynamic 

limit gives a classical radiation state in the thermodynamic limit. In 
the latter case, the result is similar but we have the parameter of the 
coherent state increasing as N, i.e. faster. Again, we get a classical 
radiation field due to the fact that the vacuum fluctuations are washed 
out in both limits. These have been amplified to the classical level and 
we have produced an intense radiation field. We have a QAMP. 

Higher order corrections can now be computed by the dual Dyson 
series as usual by the Hamiltonian 

H F = U F (t)AS z U F (t) 

that gives us 

H F = E (h,S\S x ]\AS z \[n-S,S x })\[n;S',S x ])([n-,S,S x }\ (80) 

n S,S',S X 

m,n q of q or 

The situation is more involved than the model of a single two-level 
atom given in secUJ but the approach is identical. Here, the main 
result is that, in the limit N — > 00, one gets a classical radiation field. 
This tends to become an exact result [27J; so, even if we started with 
the strong coupling expansion, we arrived to a non-perturbative result. 



(79) 



20 



9 Conclusions 



We have reviewed the strong coupling expansion as can be applied 
to time dependent problems in quantum mechanics. This approach 
proved to be very fruitful for the study of a quantum system in differ- 
ent regime, that is, in different regions of the parameter space of the 
Hamiltonian. 

Having introduced the renormalization group method for removing 
secularities in the perturbation series, in a formulation due to Kuni- 
hiro, we have built an algorithm for doing perturbation theory, making 
very simple the computation of higher order terms in the series, with- 
out any unbounded term in time. 

We have seen the method in action by the analysis of a two-level 
atom in a single radiation mode. We have obtained the Rabi oscillation 
in the strong coupling regime that have been recently observed in 
Josephson junctions. 

The generalization of this model to N two- level atoms gives a pos- 
sible description of a new effect that can be seen as a quantum ampli- 
fier (QAMP) of vacuum fluctuations of the radiation field. The effect 
appears in the thermodynamic limit N — > oo. 

We can conclude that a fruitful approach for doing perturbation 
theory is now available to analyze quantum systems, in the time do- 
main, in different regions of the parameter space of the Hamiltonian. 
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